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Abstract. In these lectures we give, first, the model-independent analysis of the ex¬ 
clusive rare decays B —> Kll and Be —> D{D*)ll with special emphasis on the cascade 
decay Be —> -D*(—> Dn)ll. We derive a four-fold angular decay distribution for this 
process in terms of helicity amplitudes including lepton mass effects. The four-fold 
angular decay distribution allows to dehne a number of physical observables which 
are amenable to measurement. Second, we calculate the relevant form factors within a 
relativistic constituent quark model, for the first time without employing the impulse 
approximation. The calculated form factors are used to evaluate differential decay rates 
and polarization observables. We present results on a set of observables with and with¬ 
out long-distance contributions, and compare them with the results of other studies. 


1 Introduction 


The flavor-changing neutral current transitions B K+X and Be D{D*)+X 
with X = 7 , l~^l~ ■, i/iy are of special interest because they proceed at the loop level 
in the Standard Model (SM) involving also the top quark. They may therefore 
be used for a determination of the Cabibbo-Kobayashi-Maskawa (CKM) ma¬ 
trix elements Vtq {q = d,s,b). The available experimental measurements of the 
branching ratio of the inclusive radiative iJ-meson decay 


Br {B ^ Xsj) = 


(3.11 ± 0.80(stat) ± 0.72(syst)) x IQ-^ 
(3.36 ± 0.53(stat) ± 0.42(syst)t!( j^(th)) 
(3.21 ± 0.43(stat) ± 0.27(syst)t;(;(^(th)) 


ALEPH 0 
X 10"'* BELLE @1 
X 10-4 CLEO II 


are consistent with the next-to-leading order prediction of the standard model 
(see, e.g. [|| and references therein): 

Br(S ^ X^7)sM = (3.35 ± 0.30) x lO’^ . (1) 


The decay B ^ K l^l {I = has been observed by the BELLE Collabora¬ 
tion 0 with a branching ratio of 


Br {B Kl+r) = [Q.lhtofi ± 0.09) x 10"® . 


( 2 ) 




2 


M.A. Ivanov and V.E. Lyubovitskij 


The recent observation of the i?c meson by the CDF Collaboration at Tevatron 
in Fermilab raises hopes that one may also explore the rare decays of the 
bottom-charm meson in the future. 

The theoretical study of the exclusive rare decays proceeds in two steps. First, 
the effective Hamiltonian for such transitions is derived by calculating the leading 
and next-to-leading loop diagrams in the SM and by using the operator product 
expansion and renormalization group techniques. The modern status of this part 
of the calculation is described in the review (and references therein). Second, 
one needs to evaluate the matrix elements of the effective Hamiltonian between 
hadronic states. This part of the calculation is model dependent since it involves 
nonperturbative QCD. There are many papers on this subject. The decay rates, 
dilepton invariant mass spectra and the forward-backforward asymmetry in the 
decays B ^ K l^l~ (Z = e, fi, t) have been investigated in the SM and its super- 
symmetric extensions by using improved form factors from light-cone QCD sum 
rules j|] . An updated analysis of these decays has been done in Q by including 
explicit 0{as) and HgcD/m.;, corrections. The invariant dilepton mass spectrum 
and the Dalitz plot for the decay B ^ K l~^l~ have been studied in by using 
quark model form factors. The B K l~^l~ decay form factors were studied 
via QCD sum rules in [p^ and within the lattice-constrained dispersion quark 
model in . Various aspects of these decays were discussed in numerous papers 
by Aliev et al. |]^. The exclusive semileptonic rare decays B Kl^l~ were 
analyzed in supersymmetric theories in m- The angular distribution and CP 
asymmetries in the decays B —> KTTe~^e~ were investigated in [^. The lepton 
polarization for the inclusive decay B —> Xsl^l~ was discussed in [|^ and | |l^ . 
The rare decays of Be D{D*)l^l~ were studied in Q by using the form 
factors evaluated in the light front and constituent quark models. 

We employ a relativistic quark model 19| to calculate the decay form 
factors. This model is based on an effective Lagrangian which describes the 
coupling of hadrons H to their constituent quarks. The coupling strength is 
determined by the compositeness condition Zh = 0 l20i|2l| where Zh is the 
wave function renormalization constant of the hadron H. One starts with an 
effective Lagrangian written down in terms of quark and hadron fields. Then, by 
using Feynman rules, the S-matrix elements describing the hadronic interactions 
are given in terms of a set of quark diagrams. In particular, the compositeness 
condition enables one to avoid a double counting of hadronic degrees of freedom. 
The approach is self-consistent and universally applicable. All calculations of 
physical observables are straightforward. The model has only a small set of 
adjustable parameters given by the values of the constituent quark masses and 
the scale parameters that define the size of the distribution of the constituent 
quarks inside a given hadron. The values of the ht parameters are within the 
window of expectations. 

The shape of the vertex functions and the quark propagators can in principle 
be found from an analysis of the Bethe-Salpeter and Dyson-Schwinger equations 
as was done e.g. in |^. In this paper, however, we choose a phenomenological 
approach where the vertex functions are modelled by a Gaussian form, the size 
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parameter of which is determined by a fit to the leptonic and radiative decays of 
the lowest lying charm and bottom mesons. For the quark propagators we use the 
local representation. In the present calculations we do not employ the so-called 
impulse approximation used previously |]^. The numerical results obtained with 
and without the impulse approximation are close to each other for light-to- 
light and heavy-to-heavy transitions but differ considerably from one another 
for heavy-to-light transitions as e.g. in the B ^ tt transitions. 

We calculate the form factors of the transition B ^ K and use them to eval¬ 
uate differential decay rates and polarization observables with and without long¬ 
distance contributions which include the lower-lying charmonium states accord¬ 
ing to [^. We extend our analysis to the exclusive rare decay Be D{D*)ll. We 
derive a four-fold angular decay distribution for the cascade B^ £>*(—*■ DTr)ll 
process in the helicity frame including lepton mass effects following the method 
outlined in Q. The four-fold angular decay distribution allows one to define a 
number of physical observables which are amenable to measurement. We com¬ 
pare our results with the ones of other studies. 


2 Effective Hamiltonian 


The starting point of the description of the rare exclusive decays is the effec¬ 
tive Hamiltonian obtained from the SM-diagrams by using the operator product 
expansion and renormalization group techniques. It allows one to separate the 
short-distance contributions and isolate them in the Wilson coefficients which 
can be studied systematically within perturbative QCD. The long-distance con¬ 
tributions are contained in the matrix elements of local operators. Contrary to 
the short-distance contributions the calculation of such matrix elements requires 
nonperturbative methods and is therefore model dependent. 

We will follow Refs.Q in writing down the analytical expressions for the 
effective Hamiltonian and paper in using the numerical values of the input 
parameters characterizing the short-distance contributions. At the quark level, 
the rare semileptonic decay b —> s{d)l~^l~ can be described in terms of the 
effective Hamiltonian: 


Hefi 


Gf 

V2 


10 


( 3 ) 


where Xt = is the product of CKM elements. For example, the standard 

set 0 of local operators for b sl^l~ transition is written as 


Ql = {SiCj)v-A {Cjbi)v-A , 

Qa = {sb)v-A J2ql9q)v-A , 

Qa = {sb)v-A Y.q{d(l)v+A , 

Ql = (1 + ^5)b, 

Qg = s^{sb)v-A{ll)v 


Q2 = {sc)v-A{cb)v-A ) 

Qa = {Sibj)v-A Y,q%<li)v-A , 

Q& ~ {Sibj)v—A ^2qidjQi)v+A , 

Qa = ^mbSia^^'^'(l + ^ 5 )T^jbj , 
Qw = ■^{sb)v-A{ll)A 


( 4 ) 
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where and F^i, are the gluon and photon field strengths, respectively; 
are the generators of the SU (3) color group; i and j denote color indices (they 
are omitted in the color-singlet currents). Labels {V ± A) stand for 7'^(1 ± 7 ®). 
( 3 i ,2 are current-current operators, Qa-e are QCD penguin operators, are 
’’magnetic penguin” operators, and Qgpo are semileptonic electroweak penguin 
operators. 

The effective Hamiltonian leads to the free quark b —> sl^l~-decay amplitude: 

Mib ^ se+n = ^ At {Cf {sb)y_^ {Il)y + Cw {sb)y_^ (II)^ (5) 

- {sia^^il + ^^)q‘'b) {Il)yY 

where Gy® = Ct — — Ce- The Wilson coefficient G|® effectively takes into 

account, first, the contributions from the four-quark operators Qi (i=l,..., 6 ) and, 
second, the nonperturbative effects coming from the cc-resonance contributions 
which are as usual parametrized by a Breit-Wigner ansatz : 


Gg® = CqFCo{ h{mc, s) + 


Stt 


E 


Vi='ip{ls),-ip(2s] 


)mvi 

mvi — imvi 


- ih(l, s) ( 4 G 3 + 4 G 4 + 3 G 5 + Ge) 

- ih(0, s) (G3 + 3G4) + ^ (3G3 + G4 + 3G5 + Ge) . 


( 6 ) 


where Gg = 3 Gi -I- Gg -I- 3G3 -I- G4 -|- 3G5 -|- Gg, ihc = rric/mB, s = q^/m?^ and 
K = I/Gq. Explicit expressions for the function hirhc, s), mi, = mb{^i) and 
can be found in Refs. 0|. The numerical values of the input parameters are 
taken from and the corresponding values of the Wilson coefficients used in 
the numerical calculations are listed in Table |^. 


3 Form factors and differential decay distributions 

We specify our choice of the momenta as pi = P 2 + ki + /c 2 with pi = ml, 
p\ = and kl = = jj? where ki and fcg are the and l~ momenta, and 

TOi, m 2 , p are the masses of initial and final mesons and lepton, respectively. 
We define dimensionless form factors by 

< K{D){p2) I s{d)-i^b\B{B^)(pi) >= F+{q^)P^+ F_{q^)q^ , (7) 

< K{D){p2)\s{d)ia^yq'' b\B{Bc){pi) >= -^ q^ Friq^), 

mi + m 2 


mi + m 2 


i < D*{p2,e2) \ d0^b\Bc{pi) >= 
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+ Pfj,Py A+{q^) P q^,P^, A_{q^) + l£^,.a/3P“/, 

i < D*{p2,e2) I dia^^q''{l + 75 ) &| Bc(pi) >= 

= gjlu Pg ao(g^) - Pit Pi' a+i<f) - i£^Li^ai3P°'q^ g(<?^)} 

where F = Pi+P 2 , q = pi-p 2 , P^; = PiJ. - qp.Pq/q^, = 9i^<'- ^ 

and el is the polarization four-vector of the D*. 


Table 1. Central values of the input parameters and the corresponding values of the 
Wilson coefficients used in the numerical calculations. 


mw 

80.41 GeV 

Gi 

-0.248 

mz 

91.1867 GeV 

C 2 

1.107 

sin^ Ow 

0.2233 

C 3 

0.011 

me 

1.4 GeV 

Cl 

-0.026 

mt 

173.8 GeV 

Cb 

0.007 

^6,pole 

4.8 GeV 

Ce 

-0.031 


mb^pole 

Cf 

-0.313 

dgcD 

0.220 GeV 

C 9 

4.344 

a~^ 

129 

Cio 

-4.669 

Osimz) 

0.119 

Co 

0.362 

\VlVtb\ 

0.0385 



\VtVtb\ 

0.008 



\VlVtb\/\Veb\ 

1 




The matrix elements of the exclusive transitions B ^ Kll and Be D{D*)ll 
are written as 

M {B{Be) - K{D*)Il) = ^ ^ {Tt M + Tt (^7^75/)} (8) 

where the quantities Tf are expressed through the form factors and the Wilson 
coefficients in the following manner: 

(a) B{Bc) K{D)ll-deca.y: 


Tt = Ptt P^ + q^ 


(* = 1,2). 


= Cf F+ + Cf Ft , 

mi + m 2 

j:iy ^ (jeff p_ _ peff ‘^rUb Pq 

mi + m 2 


ip 


(9) 
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(b) Be D*ll-deca.y: 




(* = 1 , 2 ), 


( 10 ) 


rr = - - - l-Pqg^^’^ I/W | 

TOl + TO2 1 J 




A'‘> =Cf >1„ + C5 ”o„ 


,gg 2mb{mi + m 2 ) 


= Cg^® A+ + Cf a+ 


2mb{mi + m 2 ) 


= Cf + Cf (a, - „+) ^ , 

= Cio y, = Cio ^0, = Cio . 

Let us first consider the polar angle decay distribution differential in the mo¬ 
mentum transfer squared q^. The polar angle is defined by the angle between 
q = Pi — P 2 and ki (/+Z“ rest frame) as shown in Fig. |^. One has 


(fP 

dq^d cos 9 


|P2k 

(27r)3 4mi 




gf a|^ gjp^ . 1 

(27r)3 2 7r 8m^ 2 




(11) 


- ■ (<?" l^iT + ( 9 " - V) ^ 2 ) + Lf, • (l?iT + 1^2T)} 

where |p 2 | = A^/^(mf, m^, ( 7 ^)/ 2 mi is the momentum of the final meson and 
u = y/1 -v/g' is the lepton velocity both given in the i3(i3c)-rest frame. We 
have introduced lepton and hadron tensors as 


( 12 ) 

^ij — -^i ■ 


4 Helicity amplitudes and two-fold distributions 

The Lorentz contractions in Eq. (0 can be evaluated in terms of helicity ampli¬ 
tudes as described in [^. First, we define an orthonormal and complete helicity 













Exclusive rare decays of B and Be mesons in a relativistic quark model 


7 


basis e^(m) with the three spin 1 components orthogonal to the momentum 
transfer q^, i.e. = 0 for m = ±,0, and the spin 0 (time)-component 

m = t with e^(<) = q^ j 

The orthonormality and completeness properties read 

= 9m,n , e^Lim)el{n)gjnn = g^u (13) 

with (m,n = t, ±,0) and g^n = diag( + , — , — , —). We include the time 
component polarization vector e^(t) in the set because we want to discuss lepton 
mass effects in the following. 

Using the completeness property we rewrite the contraction of the lepton and 
hadron tensors in Eq. & according to 

= (m, n)gmni'gnn'H^^ {m', n') (14) 

where we have introduced the lepton and hadron tensors in the space of the 
helicity components 


L('=)(m,n) = e'^(m)e^"(n)L^t^ im,n) = e^f^{m)e''in)H%. (15) 

The point is that the two tensors can be evaluated in two different Lorentz 
systems. The lepton tensors (m, n) will be evaluated in the ll-CM system 
whereas the hadron tensors {m, n) will be evaluated in the B{Bc) rest system. 
In the B{Bc) rest frame one has 


Pi = {rni, 0, 0, 0), 

P^ = (E2, 0, 0, -|p2|), (16) 

g'® = (go, 0, 0, IP 2 I), 

where i ?2 = (mf+m^ —g^)/2mi and go = (wi —m2+g^)/2mi. In the i3(i?c)-rest 
frame the polarization vectors of the effective current read 


(t) = Ip^I), 


e^(±) = ^(0, Tl, -i, 0), 

e^(0) = -^(|P 2 |, 0, 0, go)- 


(17) 


Using this basis one can express the components of the hadronic tensors through 
the invariant form factors defined in Eq. ( 0 )- 
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(a) B{Bc) K{D) transition: 

W^im,n) = (e'f^(m)T;) • (e^'^(n)Tjy = H\m)H^^ (n) (18) 

The helicity form factors H^{m) are given in terms of the invariant form factors. 
One has 


H\t) 

H\±) 

H\0) 


1 


{PqBl + q^r), 


= 0 , 

2mi |p2| 

■ " 


(19) 


(b) Be D* transition: 


W^{m,n) = e'<f^{m)e''in)Hli = e 1 '''(TO)e'"(n)T* - 5 “^ + 


ij — 


^a/3 


CK d '' 

P 2 P 2 

ml 


T, 


■ti 


I3u 


= (m)e'^ (r)e^ (s)^^*!)?' 






( 20 ) 


= e 


t^(m)ey(r)T;„ • (e'^''(n)ey(s)T;(^) = iJ*(m)iJt^(n). 


From angular momentum conservation one has r = m and s = n for m, n = ±, 0 
and r, s = 0 for m,n = t. For further evaluation one needs to specify the helicity 
components e 2 (m) (m = ±,0) of the polarization vector of the D*. They read 


e^(±) = 4s(0 . ±1: -fo 0 ). 4(0) = —(IP 21 , 0 , 0 , -ifa ). (21) 

V2 m 2 

They satisfy the orthonormality and completeness properties: 

(^\^{r)e2^l{s) =-5rs, e2^(r)e5^(s)5^s =(22) 

m 2 

Finally one obtains the non-zero components of the hadron tensors 

H\i) =et'^(t)ey(0)T;„ = ^ (Pg (-A), + A^) + , 

^ mi+m2m2\/q^ 

il*(±) = etM(±)ey(±)T;„ = -(-PgA(,±2mi|p2|foO> (23) 

^ mi + m 2 

H\0) = e^^(0)ey(0)T;„ 

1 1 
mi + m2 2m2i/W 


{—Pq {ml + ml 


4)^0 + 4m? |P2p^V) • 
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The lepton tensors (m, n) are evaluated in the ll-CM system ki +k 2 = 0. 
One has (see Fig. ^ 


= (V^, 0, 0, 0), 

k^ = {Ei, |ki| sin^cosy, |ki| sinffsiny , |ki|cos0), (24) 

k 2 ={Ei, — |ki| sin0cosx , — |ki| sin^siny , — |ki|cos0), 


with El = i/^j2 and |ki| = \/q^ — 4/i^/2. The longitudinal and time com¬ 
ponent polarization vectors in the II rest frame can be read off from Eq. ( p^ 
and are given by e'^(O) = (0,0, 0,1) and e{t) = (1,0, 0,0) whereas the transverse 
parts remain unchanged from Eq. (p^. 

The differential {q^,cos9) distribution finally reads 


dr{H,n ^ Hfll) _ 3 


dq'^d{cos0) 



'dEF 

2 1 

. dq^ 


, drr 

dq^ 

dq^ 




dq^ 


3.2.1 

- sm 0 ■ - 
4 2 


V ■ - cos 9 ■ , „ 

4 dq^ 

3 . 2 . 1 3 

-- sm 9 ■ - —rE -- (1 

4 2 dq^ 8 ^ 


cos 


2 2 
Integrating over cos 9 one obtains 
dr{H,„ ^ Hfll) _ 


3 , 2 . 

^ - sin^ 9 ■ 


cos^ 9) ■ 

dPr 


dq^ 


^ 1 drf 

dq^ 4 dq^ 


dq^ 








dEl^ ^ drl^ 


dq^ 


dq^ 


(25) 


(26) 


2 dq'^ 


dr^\ldrl^ 


dq^ 


2 dq^ 


drP 

dq^ 


1 drf 

2 dq^ 


where the partial helicity rates dE'^/dq^ and dE'^/dq^ {X = U,L,P,S; i,j = 
1 ,2) are defined as 


d£^ ^ Gp f a|At| \^ |p2|g^t^ 

dq'^ (27r)3 \ 27r y 12 

The relevant bilinear combinations of the helicity amplitudes are defined in Ta¬ 
ble |. 


dg2 


2 // 
^2 


dr^l 


g2 d(jf2 


(27) 
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Table 2. Bilinear combinations of the helicity amplitudes that enter in the four-fold 
decay distribution Eq. (§)• 


Definition 

Property 

Title 

Wjj = Re -b Re 

H% = Im + Im 

Hu — J^ij 

Tjij _ 

^lU ^lU 

U npolarized-transverse 

= Re - Re [hLhIJ) 

H% = Im {W+HX^) - Im 

Hy = Hy 

hYp = -Hix 

Parity-odd 

Hijl = Re {HXh]J) 

H% = Im {hXhX^) 


Transverse-interference 

HiX = Re 

H}i = Im 

Hy^ny 

Hyp = -ny^ 

Longitudinal 

= 3Re {HiHp) 

H% = 3lm {HiHl^) 

ny-ny 

HYs = -HiX 

Scalar 

HX\ = Re{mHl^) 

H%i^ = lm{HiHl^) 


Scalar-Longitudinal- 

interference 

ny = i [Re + Re 

ny = i [Im {HiHl^) + Im {hIhI^)] 


transverse-longitudinal- 
interference 

hH = i [Re {HiHl^) - Re {HIHI^)] 
H^a = 1 [Im {HXHli) - Im {hLhI^)] 


parity-Asymmetric 

= \ [Re {HXHV) + Re {H^-H^)] 

H%t = i [Im {hXhI^) + Im (Jli Jfp)] 


Scalar-Transverse- 

interference 

HiXA = I [Re {HXHl^) - Re {HtHl^)] 

HYsa = I [Im (hXhP) - Im 


Scalar-Asymmetric- 
interference 
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5 The four-fold angle distribution in the cascade decay 

Be ^ D*{^ D7r)Il. 


The lepton-hadron correlation function reveals even more structure 

when one uses the cascade decay Be —> D*{-^ Dtt)11 to analyze the polarization 
of the D*. The hadron tensor now reads 

^ (28) 

where 3°“^ (^ 2 ) = — 5 ““ +P 2 P 2 /™2 i® the standard spin 1 tensor, p 2 = ps +P 4 , 
P 3 — Pi — P 3 and p 4 are the momenta of the D and the tt, 

respectively. The relative configuration of the (H, tt)- and (/Z)-planes is shown in 

Fig. |. 

In the rest frame of the D* one has 


P 2 =(™d*,0), (29) 

P 3 = (F^n , IpsI sin 0 * , 0 , -IpsI COS 0 * ), 

= (F;^ , -IpsI sin 6 »* , 0 , IpsI cos 6 »* ), 

IpsI = A^/^(m|j.,TO|,,TO 2 )/( 2 mD.). 

Without loss of generality we set the azimuthal angle y* of the (FI, 7 r)-plane to 
zero. According to Eq. (El the rest frame polarization vectors of the D* are 
given by 

e^(±) = ^(0, ±, -z, 0), e^(0) = (0, 0, 0, -1). (30) 

The spin 1 tensor S'"" (^ 2 ) is then written as 


^““'(p2) = -5““' + ^= E e^(m)et“'(m) (31) 

ruo 

^ m =±,0 

Following basically the same trick as in Eq. (^^ the contraction of the lepton 
and hadron tensors may be written through helicity components as 




fl' V' 




(32) 


L^{m,n)gmm'gnn' (e'''''(TO')eJ“(r)r^„) {n )el°‘{s)Tl^ 
^ 3Br(F)* ^ Dtt) 

P3(^2{r) ■ p^e'^is) -——- 

2 iPsI 


3Br(£>* ^ Dtt) 
2|p3| 




^ L'"{rn,n)W(rn)H'^^(n) ■ P3e2{m) ■ pselin) 

m , n =±,0 
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- L^(t,n)W{t)H'^\n) ■ P3e2i0) ■ pselin) 

n—±,0 

- Y ■ P3e2im) ■ pseliO)] 

m =±,0 

Using these results one obtains the full four-fold angular decay distribution 
dq^dcosed{x/2Tr)dcose* ^ ^ ^ ^ 


drh\ 

dr^^ 

dq^ 

dq^ 

drr\ 


dq^ ) 



drp^ 

dq^ 

dq^ 


3.2., a 3.2.,*! dB 

-- sin^ 6» • cos 2x • 7 sin^ d ■ - [ —^ -|- —^ 
4 ^4 2\dq^dq 


9 • o.. 1 fdr}^ drp 

+„ [-2cos9- 

4 4 dq^ 

9 . . -on* 1 drf\ 

--smd.cosx.sm20 • - ^j 

9 * fdrjf drf}\ 

+ -sind.sinx.sin2d + 

9 . o. .. o.* fdrii drn\ 

--sm2d.smx .81020 . ^ j 


3 9 3 9 ,1 dPrY 3 , 9.3 2 , df.f 

-1-7 sin^ 0.7 sin^ • 7 ' -7^ “ E (^ + <^0® d) ■ - sin^ 0 . -7^ 

4 4 2 dq^ 8^ ^4 dq^ 

3 2 ., 3 2 .,* 1 3 . 2 ., 3 2 .,* drl^ 

-\— cos^ 0 . - cos^ 0 . - . —77-sm^ 0 . - cos^ 0 . , H, 

2 2 2 dg 2 4 2 dg 2 


3 3 

-b- sin^ 0 . cos 2x . - sin^ 9* 

4 ^4 


dr“ dPp 
~df ^ ~df 


9 . on . on* 3 2n* ^ dPf 

-- sm20 . cosx • Sm20 • 7 — 77 -+- 7 - 7 - + 7 cos^ 0 . 7 - 74 - 

8 2 \ dq^ dq^ I 2 4 dq^ 


-\ -sin 20 . sin x . sin 20* . 

16 

9 

-sin^ 0 . sin 2x . sin^ 0* 

16 


.7711 

IA I IA 

dq^ dq^ 2 

(YY.'Ejt 

\ dq^ dq^ 
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Integrating Eq. over cos 0* and x one recovers the two-fold , cos 0) dis¬ 
tribution of Eq. ( p5|). Note that a similar four-fold distribution has also been 
obtained in Refs.(|14[,p5t,]2^, using, however, the zero lepton mass 

approximation. If there are sufficient data one can attempt to fit them to the 
full four-fold decay distribution and thereby extract the values of the coefficient 
functions dFx/dq^ and, in the case I = t the coefficient functions dFx/dq^. In¬ 
stead of considering the full four-fold decay distribution one can analyze single 
angle distributions by integrating out two of the remaining angles as done in 
Ref. lai. 


6 Model form factors 


We will employ the relativistic constituent quark model |^,|^ to calculate the 
form factors relevant to the decays B Kll and Be D{D*)ll. This model 
is based on an effective interaction Lagrangian which describes the coupling 
between hadrons and their constituent quarks. 

For example, the coupling of the meson F[ to its constituent quarks qi and 
q 2 is given by the Lagrangian 

Cint{x) = gHH{x) JdxiJdx2FH{x,xi,X2)q{xi)FHXHq{x2) ■ (34) 

Here, Xh and Fh are Gell-Mann and Dirac matrices which entail the flavor and 
spin quantum numbers of the meson held F[{x). The function Fh is related 
to the scalar part of the Bethe-Salpeter amplitude and characterizes the hnite 
size of the meson. The function Fh must be invariant under the translation 
Fh{x -I- a, xi -I- a, X 2 + a) = Fh{x, xi,X 2 )- 

In our previous papers we have used the so-called impulse approximation for 
the evaluation of the Feynman diagrams. In the impulse approximation one omits 
a possible dependence of the vertex functions on external momenta. The impulse 
approximation therefore entails a certain dependence on how loop momenta are 
routed through the diagram at hand. This problem no longer exists in the present 
full treatment where the impulse approximation is no longer used. In the present 
calculation we will use a particular form of the vertex function given by 

Fh{x, xi,X 2 ) = s(x- -X 2 )^)- (35) 

\ mi +m2 J 

where m-i and m 2 are the constituent quark masses. The vertex function Fh 
evidently satishes the above translational invariance condition. As mentioned 
before we no longer use the impulse approximation in the present calculation. 

The coupling constants gn in Eq. (^J) are determined by the so called com¬ 
positeness condition proposed in [^| and extensively used in |^. The compos¬ 
iteness condition means that the renormalization constant of the meson field is 
set equal to zero 
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where is the derivative of the meson mass operator. For the pseudoscalar 
and vector mesons treated in this paper one has 


n'pip") = 


'd'^k 


2p^ dp°‘ J 

X tr 




n'vip^) = I 

X tr 


7 ^ S ' i (/^ + W21 - W12 

1 . d fd^k 


r 


pPp’' 


p 


'2 2p2^“dp« J 


YSi{ji + W21 ■^)-f'^S2{^- W12 


where Wij = rrij /(m^ + rrij). 

The leptonic decay constant fp is calculated from 

+ W21 f)l^S2(}t - W12 fj 


= fpP^ 


(37) 


The transition form factors P{pi) P{p 2 ),V{p 2 ) can be calculated from the 
Feynman integral corresponding to the diagram of Fig. 0: 


A^'^{Pl,P2) = ^^^1^2^^^ f^;^3pi-ik + Wi3Pif)3p,(v){-ik + W23P2f) 


X tr 


^2(^+ ^2)T^^i(^+ fil)-f^S3{li)rout 


(38) 


where r>^ = or and Pp'y = y®, yi/C^- 

We use the local quark propagators 


S^{^) = 


1 


(39) 


TTLi— ji ’ 

where is the constituent quark mass. We do not introduce a new notation for 
constituent quark masses in order to distinguish them from the current quark 
masses used in the effective Hamiltonian and Wilson coefficients as described in 
Sec. II because it should always be clear from the context which set of masses is 
being referred to. As discussed in [|^,^, we assume that 

mp < mi + m 2 (40) 

in order to avoid the appearance of imaginary parts in the physical amplitudes. 
The fit values for the constituent quark masses are taken from our papers 
Hl^ and are given in Eq. 


ruu 


ms 


me 


mb 


0.235 0.333 1.67 5.06 GeV 


(41) 
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It is readily seen that the constraint Eq. (^) holds true for the low-lying fla¬ 
vored pseudoscalar mesons but is no longer true for the vector mesons. In the 
case of the heavy mesons D* and B* we will employ identical masses for the 
vector mesons and the pseudoscalar mesons for the calculation of matrix ele¬ 
ments in Eqs. (^^,(|3^) and It is a quite reliable approximation because of 
{mo* —rnE))/mD ~ 7% and {uib* —mB)/mB ~ !%• In this vein, our model was 
successfully developed for the study of light hadrons (e.g., pion, kaon, baryon 
octet, zl-resonance), heavy-light hadrons (e.g., D, Dg, B and B^j-mesons, Aq, 
Eq, Sq and I^g-baryons) and double heavy hadrons (e.g, J/'P, T and Bc-mesons, 
Eqq and 17 qq baryons) [|l8|Jl9|| . To extend our approach to other hadrons we 
had to introduce extra model parameters or do some approximations, like, e.g., 
to introduce the cutoff parameter for external hadron momenta to guarantee 
the fulfillment of the mentioned above ’’threshold inequality”. Therefore, at the 
present stage we can not apply our approach for the study of rare decays involv¬ 
ing K* mesons. Probably, it will be a subject of our future investigations. 

We employ a Gaussian for the vertex function A'j^) = exp(— kl/Al) 

where is the Euclidean momentum and determine the size parameters Ajj by 
a fit to the experimental data, when available, or to lattice simulations for the 
leptonic decay constants. The quality of the fit can be seen from Table The 
branching ratios of the semileptonic decays are shown in Table ^ The numerical 
values for Ah are = 1 GeV, Ak = 1.6 GeV, Ab = 2 GeV and yls = 2.25 GeV 
for all K, D and B partners, respectively. 


Table 3. Leptonic decay constants fn (MeV) used in the least-square fit. The values 
are taken either from PDG [s^ ] or from the Lattice quenched (upper line) 
and unquenched (lower line). 


Meson 

This model 

Expt/Lattice 


131 

130.7 ±0.1 ±0.36 

K+ 

161 

159.8 ± 1.4 ±0.44 

D+ 

211 

203± 14 

226± 15 

D+ 

222 

230± 14 

250± 30 

B+ 

180 

173± 23 

198± 30 

B° 

196 

200± 20 

230± 30 

B+ 

398 
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Table 4. Semileptonic decay branching ratios. 


Meson 

This model 

Expt. 

TT^ ^ 

1.03 ■ 10"® 

(1.025 ± 0.034) • 10“® 

K+ TT°l+iy 

4.62 • 10"^ 

(4.82 db 0.06) • 10"^ 

B+ D°l+u 

2.40 • 10"^ 

(2.15 ±0.22) • 10"^ 

B+ D*°l+v 

5.60 • 10"^ 

(5.3 ±0.8) • 10“^ 

B+ D°l+v 

2.05 • 10"® 


B+ ^ D*°l+v 

3.60 • 10“® 



We are now in a position to present our results for the B{Bc) K{D^ D*) 
form factors. We have used the technique outlined in our previous papers il8|Jl9| 
for the numerical evaluation of the Feynman integrals in Eq. (^). The results 
of our numerical calculations are well represented by the parametrization 


F{s) 


m 

1 — as + bs'^ ' 


(42) 


Using such a parametrization facilitates further integrations. The values of E(0), 
a and b are listed in Tables |[ 


Table 5. Parameter values for the approximated form factors 
F{s) = f’(0)/(l — as + bs^) (s = g^/m|). 


B Kll 

F— Ft 

F{0) 

0.357 -0.275 0.337 

a 

1.011 1.050 1.031 

b 

0.042 0.067 0.051 


Be D{D*)ll 

F+ 

F- 

Ft 

^0 

71+ 

A- 

V 

ao 

a+ 

9 

E(0) 

0.186 

-0.190 0.275 

0.279 0.156 

-0.321 0.290 

0.178 0.178 0.179 

a 

2.48 

2.44 

2.40 

1.30 

2.16 

2.41 

2.40 

1.21 

2.14 

2.51 

b 

1.62 

1.54 

1.49 

0.149 

1.15 

1.51 

1.49 

0.125 

1.14 

1.67 


At the end of this section we would like to discuss the impulse approxima¬ 
tion used in our previous papers [|l8yi9|] . It was simply assumed that the vertex 
functions depend only on the loop momentum flowing through the vertex. The 
explicit translational invariant vertex function in Eq. (|3^ allows one to check the 




















Exclusive rare decays of B and Be mesons in a relativistic quark model 


17 


reliability of this approximation. We found that the results obtained with and 
without the impulse approximation are rather close to each other except for the 
heavy-to-light form factors. We consider the B 7r-transition as an example to 
illustrate this point. The calculated values of the form factor at = 0 

are 

{ 0.27 exact 

0.48 impulse approximation 

One can see that the value of the form factor at = 0 calculated without the 
impulse approximation is considerably smaller than when calculated with the 
impulse approximation. Its value is close to the value of QCD SR estimates, see, 
for example, [^: F_[^’^(0) = 0.30. 

7 Numerical results 

We list our numerical results for the branching ratios in Table pj. When comparing 
the values of the branching ratios with those obtained in and 0 one finds 
that they almost agree with each other. 


Table 6 . Decay branching ratios without (with) long distance contributions. 


Ref. 

Bt(B ^ Kfi+fi-) 

Br(B ^ Kt+t-) 

Br(R ^ Kw) 



0.57- lO"*’ 

1.3 ■ 10“' 




(0.35 ±0.12) • 10^'’ 



JL 

11 

0.44 • lO"” 

1.0 ■ 10“' 

5.6- 10““ 


2 

0.5 • lO"*’ 

1.3 ■ 10“" 


our 

0.55 (0.51) • 10~“ 

1.01 (0.87) • 10~ ' 

4.19 • 10"'’ 



our 

|17| 

Br{Bc ^ ) 

Br(Bc ^ Dd 

Br(Bc ^ Ds 

Br(Bc ^ Dl 

0.44 (0.38) • 10"'' 
0.71 (0.58) • 10"® 
0.97 (0.86) • 10"^ 
1.76(1.41) • 10"'^ 

0.41 (0.33) • 10"® 
1.01 (0.78) • 10"® 
1.36(1.12) • 10"’^ 
4.09 (3.14) • lO"’’ 

Br(Bc ^ DdT^T~) 
Br{Be^D*dT+T-) 
Br(Bc ^ Ds r+r ) 
Br{Be^ DIt+t-) 

0.11 (0.09) • 10’® 
0.11 (0.08) • 10"® 
0.22(0.18) • 10"^ 
0.22 (0.15) • lO”’’ 

0.13 (0.11) • 10”® 
0.18 (0.13) • 10"® 
0.34(0.27) • 10"’^ 
0.51 (0.34) • 10"’^ 

Br(Bc ^ Dd vu) 
Br(i3c ^ Dd vv) 
Br(Bc ^ Ds vv) 
Br(i3c —> D* vv) 

3.28 • 10"® 
5.78 • 10"® 
0.73 • 10“® 
1.42 • 10~® 
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Fig. 1. Definition of angles 9, 6* and y in the cascade decay Be D*(—> Dn)U. 


Pi 


B{D,) 



P2 


K{D,D*) 


Fig. 2. Feynman diagram describing the form factors of the decay -B(-Bc) ^ K{D, D*)ll 
















